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We report inelastic time-of-flight and triple-axis neutron scattering measurements of the excitation 
spectrum of the coupled antiferromagnetic spin-1 Heisenberg chain system CsNiCla. Measurements 
over a wide range of wave-vector transfers along the chain confirm that above Tm CsNiCla is in a 
quantum-disordered phase with an energy gap in the excitation spectrum. The spin correlations fall 
off exponentially with increasing distance with a correlation length ^ = 4.0(2) sites at T = 6.2 K. 
This is shorter than the correlation length for an antiferromagnetic spin-1 Heisenberg chain at this 
temperature, suggesting that the correlations perpendicular to the chain direction and associated 
with the interchain coupling lower the single-chain correlation length. A multi-particle continuum is 
observed in the quantum-disordered phase in the region in reciprocal space where antiferromagnetic 
fluctuations are strongest, extending in energy up to twice the maximum of the dispersion of the 
well-defined triplet excitations. We show that the continuum satisfies the Hohenberg-Brinkman sum 
rule. The dependence of the multi-particle continuum on the chain wave- vector resembles that of the 
two-spinon continuum in antiferromagnetic spin-1/2 Heisenberg chains. This suggests the presence 
of spin-1/2 degrees of freedom in CsNiCls for T < 12 K, possibly caused by multiply-frustrated 
interchain interactions. 



PACS numbers; 75.25.-l-z, 75.10.Jm, 75.40.Gb 
I. INTRODUCTION 

Nearly two decades after Haldane's conjecture! the 
interest in quantum spin dynamics of low-dimensional 
antiferromagnets continues unabated. Much of the re- 
search is focused on one-dimensional (ID) antiferromag- 
nets with small spin S, where quantum fluctuations are 
strong because the classical Neel ground state is not an 
eigenstate of the Hamiltonian. The ground state of such 
an antiferromagnetic (AF) chain is a spin singlet and 
long-range magnetic order is absent, even at T = K. 
The spin correlation function in the ground state of AF 
Heisenberg chains depends on the spin quantum num- 
ber 5": for 5 = 1/2 the spin-spin correlations fall off as 
a power law with increasing distance (quasi long-range 
ordered ground state) but for S = I the spin-spin coi^ 
relations fall off exponentially with increasing distance.El 
The excitation spectrum is also different: for S = 1/2 
the elementary excitations are spin-1/2 particles called 
spinous and the spectrum extends to zero energy trans- 
fer while for S = 1 the elementary excitations are spin-1 
particles with a gap (Haldane excitations). 

Inelastic neutron scattering is one of the most direct 
methods to measure the magnetic excitations in solids. In 
neutron scattering experiments, the spin of the neutron 
changes by AS'z = 0, ±1 along a quantization axis so that 
at low temperatures the response is dominated by broad 
continua scattering if the elementary particles carry spin- 



1/2. Sharp modes are observed, however, if the elemen- 
tary particles carry spin-1. The spectra observed with 
neutron scattering are therefore fundamentally different 
in AF spin-1/2 and spin-1 Heisenberg chains: for a spin- 
1 chain at low temperatures the excitation spectrum is 
dominated by well-defined modescTLl while in the case 
of a spin-1/2 chain the dominant contribution is a con- 
tinuum consisting of two-spinon states. The two-spinon 
continuum has a characteristic wave-vector dependence 
reflecting the momentum and energy addition of two par- 
ticles and it has an upper energy boundary at the AF 
point which is double the maximum single-particle en- 
ergy and decreases away from it. The obscfyation of such 
a scattering continuum in a spin systemcl thus clearly 
suggests that spin-1/2 particles are its elementary exci- 
tations rather than spin-1 particles. 

Considering the fundamental difference between AF 
spin-1/2 and spin-1 chains, it came as a surprise when 
we observed in the spin-1 chain compound CsNiCla a 
multi-particle scattering continuum that resembles the 
two-spinon continuum of AF spin-1/2 chains.Q Its wave- 
vector dependence and its intensity suggest that it does 
not originate from two- or three-particle excitations of 
the elementary spin-1 excitations but instead from spin- 
1/2 degrees of freedom in the system. Their detailed 
study promises to shed light on the transition from spin- 
1/2 particles in quantum critical AF spin-1/2 chains and 
the spin-1 particles found in spin liquids like the gapped 
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Haldane quantum antiferromagnets- 

The focus of our previous papeiQ was on the excita- 
tion spectrum of the quantum-disordered phase at the AF 
point Qc = 1 and it was shown with three independent 
neutron scattering measurements that there is a multi- 
particle scattering continuum extending from the gapped 
onset of the excitations up to 12 meV. The present pa- 
per contains a study of the excitation spectrum for a 
wide range of wave- vector transfers and for both the well- 
defined excitations and the multi-particle continuum. 
For the well-defined excitations the results largely con- 
firm previojis results in the quantum-disordered phase of 
CsNiCla.B'Llfl Data measured with a sample mount which 
consists exclusively of aluminum to avoid any background 
scattering originating from hydrogenous materials con- 
firm the previous estimate that 12(2)% of the scattering 
at the AF point Qc = 1 is carried by the multi-particle 
continuum. A detailed description of the multi-particle 
states is given for a wide range of wave-vector transfers. 
The integrated intensity of the scattering is shown to be 
consistent with the first-moment sum rule. The charac- 
teristics of the continuum are compared to those of the 
two-spinon continuum in AF spin-1/2 Heisenberg chains. 



II. PROPERTIES OF CSNICL3 

CsNiCls is one of the most studied AF coupled spin-1 
chain systems. It crystallizes in a hexagonal crystal struc- 
ture, Dq^ space group, and the lattice constants at low 
temperatures are a — 7.14 A and c — 5.90 A (Fig. |^). 
The magnetic moments are carried by Ni^^-ions which 
interact via a superexchange interaction involving Cl^- 
ions. Because the superexchange path along the c-axis 
contains only one Cl^-ion and perpendicular to the c-axis 
two Cl^-ions (Figs. |l|, |2|), the superexchange between the 
Ni^+-ions along the c-axis is much stronger than perpen- 
dicular to the c-axis. The spin Hamiltonian of CsNiCla is 
that of a system of coupled spin-1 chains with a strong in- 
trachain interaction J and a weak interchain interaction 
J'. It can be written as 
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h^jY^s,- S,+i + j' • - ^ ^(^f )' • (1) 

i <i-j> i 

The superexchange interaction J along the c-axis was 
measured with a high-field magnetization measurement 
of the magnetic saturation field and is 2.28 meVfl The 
superexchange interaction perpendicular to the c-axis is 
J' = 0.044 meV as determined from the measurement 
of the spin-wave energies in the antiferromagnetically or- 
dered structure = 2 K and comparing them with 
spin-wave theoryBu The weak easy-axis Ising anisotropy 
D — A fJ,eY, which was derived from the spin- flop transi- 
tion field extrapolated to T = K, is small enough that 
CsNiCla is a gop4 example of an isotropic Heisenberg 
antiferromagnet .113 



Due to the interchain interaction, the coupled spin 
chain system undergoes three-dimensional (3D) long- 
range ordering below T/v = 4.84 K, where the c-axis com- 
ponents of the magnetic moments order.E2l At Tcantcd^ 
4.4 K, the components perpendicular to the c-axis order 
too. In the resulting magnetic structure, | of the spins 
are canted awa|5z_from the c-axis by 44° at 4.4 K and 
59° at 1.6 K,Ma due to the easy-axis anisotropy and 
the magnetic frustration of the triangular lattice. Ex- 
trapolated to r = K, the ordered magnetic moment is. 
1.05/iB, considerably less than the free-ion value 2(13^3 
This reduction can be ascribed to strong quantum fluc- 
tuations in the 3D ordered phase. Strong quantum fluc- 
tuations in the 3D ordered phase are also apparent in 
the excitation spectrum, where excitations derived from 
longitudinal Haldane modes were observed which_cannot 
be described by conventional spin-wave theory.Eil Their 
existence is related to longitudinal fluctuations of the or- 
dered moment, as shown in a calculation of the excita- 
tions which includes the gapped Haldane mode_£Oupled 
to transverse gapless spin-wave-like excitations.E£l 




FIG. 1. Crystal structure of CsNiCls. The Ni atoms are 
shown as dark spheres, and the Cs and CI are shown as light 
spheres. The figure shows four Ni chains on a triangular lat- 
tice running along the c-axis. The Ni atoms are coupled 
along the c-axis by an AF superexchange interaction along 
the Ni-Cl-Ni paths. 

In the quantum-disordered phase above T/v, the ex- 
citation saectrum of CsNiCla has an energy gap (Hal- 
dane gap)cl with a minimum at the AF wave-vector. The 
gap has roughly the energy expected of a Haldane gap, 
which is a consequence not of a single ion anisotropy 
but of strong quantum fluctuations. Tha. fixcitations 
are a triplet of well-defined 5" = 1 modesj3'E3 and the 
spin-spin static correlation function falls off exponentially 
with increasing distance. For wave-vectors near the nu- 
clear zone center, a broadening of the excitation spectrum 
may arise fcom the theoretically predicted two-particle 
continuum.EZI 

The magnetic susceptibility and the magnpt|ic .specific 
heat have a broad maximum near T = 30 kJIjO, which 
had been traditionally interpreted as the temperature 
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where short-range correlations set in. A neutron scat- 
tering experiment showed that the excitations at. the . AF 
point persist as a resonance up to at least 70 so 
that short-range correlations survive even to tempera- 
tures of the order of the spin band width. 




FIG. 2. Crystal structure of CsNiCls projected onto the 
basal plane. The Ni atoms are shown as dark spheres, and 
the Cs and CI are shown as light spheres. Ni atoms belong- 
ing to different chains are coupled by an AF superexchange 
interaction along the Ni-Cl-Cl-Ni paths. 



III. EXPERIMENTAL DETAILS 

An inelastic neutron scattering experiment measures 
the response of a coupled spin system as a function of en- 
ergy and momentum transfer. The magnetic scattering 
cross section is directly related to the dynamic structure 
factors S°''^{Q,uj), which arc the Fourier transform of 
the space- and time-dependent spin-pair correlation func- 
tions. Defining the neutron energy and momentum trans- 
fers to the spin system as uj — Ei — Ef and Q = ki — kf , 
respectively, the magnetic scattering cross section for un- 
polarized neutrons can be written as 



dQcLo 
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QMp)S"''iQ,u;), (2) 



where a, P — x,y, z and /(Q) is the magnetic form fac- 
tor. For a spin system with Heisenberg interactions in 
the paramagnetic phase, the entire magnetic response 
is inelastic and the dynamic structure factors have the 
property S{Q,uj) = S^^— S^^ — S^^ because the system 
has unbroken rotational symmetry. For the Hamiltonian 
of CsNiCls, the total z-component of spin is a constant of 
motion if D J2ii^i)'^ neglected {S^ — ^ ■ commutes 
with the Hamiltonian) and S"^ = for a ^ p. 

A single crystal of CsNiCla 5 x 5 x 20 mm'^ was 
used for the experiments. The experiments were per- 
formed using the triple-axis spectrometer DUALSPEC 
at the Chalk River Laboratories, Canada, and the chop- 
per time-of-flight spectrometer MARI at the pulsed spal- 
lation source ISIS of the Rutherford Appleton Labora- 
tory, United Kingdom. The sample was kept sealed in 



an aluminum can containing helium gas to prevent the 
absorption of water. For the DUALSPEC experiment, 
the sample was wired to an aluminum plate and placed 
inside an aluminum can such that no sealant could cause 
extraneous scattering of neutrons. 

The measurements with the DUALSPEC triple-axis 
spectrometer were performed at 6 K. The sample was 
mounted in a cryostat with its (hhl) crystallographic 
plane in the horizontal scattering plane. The scattered 
neutron energy was set to 14.51 meV. A graphite filter 
was used to absorb the higher order reflections from a ver- 
tically focusing pyrolytic graphite monochromator and a 
flat analyzer. The effective collimation from reactor to 
detector was 39'-48'-51'-72' and gave an energy resolu- 
tion of 0.8 meV, as determined from the quasi-elastic 
incoherent peak. The longitudinal wave-vector width of 
the (002)-peak and (220)-peak were 0.022 and 0.019 (fuU 
width at half maximum (FWHM) in reciprocal lattice 
units), respectively, and the calculated vertical resolution 
was 0.29 A-i. 
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FIC. 3. Schematic drawing of the chopper time-of-flight 
spectrometer MARI at ISIS. A monochromatic neutron beam 
is generated using a chopper. The neutrons are scattered at 
the sample and then detected in the detector banks located 
at a distance of 4.02 m in a half circle below the sample. 
Depending on the incident energy and the orientation of the 
sample, different wave-vector and energy transfer regions are 
accessed during an experiment. 

The time-of-flight spectrometer MARI produced a 
monochromatic beam by means of the Fermi chopper 
spinning at 150 Hz. The scattered neutrons were de- 
tected by an array of He'^ tube detectors typically 30 cm 
long (perpendicular to the principal scattering plane) and 
with a diameter of 2.5 cm situated 4.02 m from the sam- 
ple position. The low-angle bank consists of 136 detec- 
tors arranged symmetrically around the incident neutron 
beam direction and having scattering angles from 3.43 to 
13.29°. The detectors in the high-angle banks are ar- 
ranged in three strips vertically below the sample and 
cover the scattering angles from 12.43 to 134.14°, where 
the detector tubes are aligned perpendicular to the ver- 
tical scattering plane. The central strip is in the verti- 
cal scattering plane, and the two side strips of detectors 
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are out of the vertical plane by an average of 5.56° for 
high scattering angles. At small scattering angles, the 
detectors in the two side strips are out of the principal 
scattering plane by up to 21.36° and the neutrons de- 
tected in these detectors have a considerable wave- vector 
component perpendicular to the scattering plane. 

The experiment using MARI was performed with the 
sample oriented such that the (hhl) crystallographic 
plane was in the vertical scattering plane and the c-axis 
was perpendicular to the incident beam direction. The 
incident neutron energy Ei was set either to 20 or 30 meV 
and the energy resolution was 0.35 and 0.4 meV, respec- 
tively, as determined from FWHM of the quasi-elastic in- 
coherent scattering peak. The resolution in wave-vector 
transfer at zero energy transfer was typically 0.02 A^^ 
along the c*-axis and along the [110] direction and up 
to 0.19A^^ perpendicular to the scattering plane if only 
the central detector bank was used. Both the energy and 
the wave- vector resolution improved with increasing en- 
ergy transfer. The measurements were performed with 
the sample at temperatures between 6.2 and 12 K. The 
scattering was measured for a total proton charge be- 
tween 4000 and 8600 /iAhr at an average proton current 
of 170 /lA. 

The simultaneous measurement of continuous energy 
transfers lu together with the large range of scattering 
angles allows the measurement of the dynamical struc- 
ture factor S{Cl,Lu) over the scattering surface, which is 
a two-dimensional subspace of (Qa, Qa w), where both 
Qa and Qc lie in the scattering plane. The scattering sur- 
face may then be projected down onto {Qc, a;)-plane and 
scans with a constant energy transfer or a wave-vector 
transfer along the chain can be obtained by binning the 
data appropriately. These scans will henceforth be called 
constant-w and constant-Qc scans, respectively. In our 
case, because CsNiCls is a ID magnet and the interchain 
dispersion is sufficiently small, S{Qc,uj) mainly reflects 
the excitation spectrum of the spin chains. 

IV. EXPERIMENTAL RESULTS 
A. Haldane excitations 

The quantum-disordered magnetic phase of CsNiCla 
at temperatures between 6.2 and 12 K was studied for 
a wide wave-vector and energy range using the MARI 
time-of-flight spectrometer. Fig. ^ shows three constant- 
Qc scans at 6.2 K obtained using the binning procedure 
described above. Only data measured in the central de- 
tector bank are shown in the Fig. |4[ This is because 
at low scattering angles the two side detector banks are, 
as mentioned, considerably out of the vertical scatter- 
ing plane and so the dynamic structure factor S{Q, lu) is 
different for the three detector banks. 

The data are treated as if there was a well-defined ex- 
citation and some extra scatteringcl but the data could 



also be interpreted as a continuum with a well-defined 
onset at the low-energy boundary. The well-defined ex- 
citation has the most intensity at the AF point Qc = 1, 
which corresponds to a tt phase difference between the Ni 
spins along the chains. The intensity decreases towards 
Qc = 0.5, but it is visible down to Qc = 0.2. The scatter- 
ing was fitted with an antisymmetrized Lorentzian form 
weighted by the Bose factor 

S{Q,uj) = A-{n{uj) + l)- 

((cc;-e(Q))2+r2 ~ (co + e(g))2+r2) ' 

and convoluted with the line-shape of the quasi-elastic 
incoherent scattering, which is a good estimate of the res- 
olution particularly at low energy transfers. Here, e{Q) is 
the dispersion relation, F is the excitation width, n(Lu) + l 
is the population factor, uj is the neutron energy trans- 
fer and A is a scaling constant. For the excitations close 
to Qc = 1, we found that a Lorentzian line-shape gave 
fits with consistently lower sum of discrepancies than 
for an antisymmetrized Gaussian cross-section, in agree- 
ment with opt [Measurements obtained using the RITA 
spectrometeind'Ell. 
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FIG. 4. Neutron scattering intensity at 6.2 K as a function 
of energy transfer for three different wave- vector transfers Qc 
along the chain. The data were measured using the central 
detector bank of the MARI time-of-flight spectrometer and 
the wave-vector range sampled for each Qc in ±0.025. The 
incident energy Ei was 30 meV. The solid line is a fit as 
described in the text and the dotted line shows the onset of 
the two-magnon continuum as predicted by the NLcrM. 
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Fits were performed to obtain the intensity, the fre- 
quency and the width of the well-defined excitations. 
An energy-independent background increasing with the 
wave- vector transfer Qc along the chain was chosen such 
that the scattering at high energy transfers matches the 
scattering at low energies below the well-defined exci- 
tation. The excitation energy and width are shown in 
Fig. I The data for Qc > 1.6 were not analyzed due to 
their contamination with phonons. The maximum of the 
dispersion is 6.19(9) meV, in agreement with previous 
measurementffl and the excitation energy at = 1 is 
1.59(5) meV. 
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FIG. 5. Excitation energy as a function of wave-vector 
transfer along the chain. The solid line is a fit to the exci- 
tation energy as described in the text. Left inset: Lorentzian 
half width F as a function of wave-vector transfer Qc along 
the chain. Right inset: Wave- vector transfer along [Qa, Qa, 0] 
and Qc calculated from the fitted energy transfer and Qc- 



The observed dispersion of a sing, 
often written phenomcnologically as 
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oj{Qc) = ^/A2+t;2sin^(Qe7r)+a2 cos2(^7r), 



(4) 



where A is the Haldane gap, is the spin velocity and a 
determines the asymmetry of the dispersion about Qc = 
0.5. For a two-particle gap of 2 A as Qc ^ one would 
expect a — VSA. Because the chains are coupled in 
CsNiCla the excitations have a dispersion perpendicular 
to the chain direction. The coupling of the chains can be 
taken into account in a Random Phase Approximaiipn 
(RPA) and the dispersion relation can be written asE3 

c^2(Q) ^ ^(gj . {uj{Qc) + 2 • J'iQa, Qb) ■ S{Q)) . (5) 

S{Q) is the static structure factor and J'{Qa,Qb) is the 
Fourier transform of the exchange interaction perpendic- 
ular to the chain axis. The Fourier transform J'{Qa, Qa) 
along the reciprocal [1, 1, 0]-direction, which is of rele- 
vance for the present experiment, is given by 



J'iQa, Qa) = J' • (4 cos(27rQa) + 2 cos(47rQ„)) , (6) 

and vanishes for Qa = 0.19. The static structure factor 
falls rapidly as jQc — 1| increases, weakening the effect of 
the interchain exchange interaction on the dispersion. 

The solid line in Fig. || corresponds to a fit of Eq. ^ 
to the observed dispersion taking into account the wave- 
vector dependence of the measured excitation perpendic- 
ular to the chain (inset of Fig. ||) and fixing the excitation 
band width pcrpcndiciilar to the chain axis by setting 
a;(|, |, 1) = 0.4 meV.Ll The fitted parameters are A = 
1.24(4) meV 0.54(2)J, i;^ = 5.70(7) meV = 2.50(3)J 
and a = 2.5(3) meV = 1.1(1)J = 2.0A. ^ 

A is consistent with our previous measurementsEj and 
is higher than the gap enesgir of 1.14 meV calculated for 
the non-linear sigma modelcJ (NLcrM) which includes the 
temperature renormalization at 6.2 K. This is possibly 
due to an additional upward renormalization of tiie gap 
energy caused by the coupling of the spin chains .t3 The 
observed energy a.t Qc = I, OJ = 1.59(5) meV, is higher 
than A because it corresponds to an excitation with a 
different wave- vector transfer perpendicular to the chain 
axis. For weakly-coupled spin chains A can be directly 
measured at Q = (0.19, 0.19, 1) where the Fourier trans- 
form of the interchain coupling vanishes. In this experi- 
ment, however, the gap excitation could only be sampled 
at Q = (0.16, 0.16, 1) as shown in the right inset of Fig.||, 
and here the excitation energy is higher. 



3 

CO 



03 



AF spin-1 chain is ^ 



• Q„=0.25 (r.Lu 
o Q„=0.75 (r.Lu 



CsNiClj T=6.2K 




4 5 6 

Energy (meV) 

FIG. 6. Neutron scattering intensity as a function of energy 
transfer for Qc ~ 0.25 and Qc = 0.75. The energy of the exci- 
tation at Qc = 0.25 is clearly higher than at Qc = 0.75. Inset: 
The solid circles represent the excitation energy for Qc > 0.5 
(bottom axis) and the open circles data for Qc < 0.5 (top 
axis), both corrected for the Qa-dependence of the dispersion. 
Hence the dispersion is not periodic about Qc ~ 0.5. 

The spin velocity Vs is in excellent . agre ement with nu- 
merical results predicting Vg = 2.5 J.cjEh a is consistent 
with previous results obtained for NENFq and CsNiClsEEl 
and its non-zero value reflects the asymmetry of the dis- 
persion with respect to reflection at Qc — 0.5 as would 
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occur for a simple antiferromagnet. Fig. ^ shows the 
asymmetry about Qc = 0.5 as measured using the MARI 
spectrometer. The figure compares two equivalent scans 
at Qc ~ 0.25 and Qc — 0.75, and the excitation at 
Qc = 0.75 clearly has a lower energy than the one at 
Qc ~ 0.25. Due to the dependence of the excitations on 
Qa (right inset of Fig. H), the energy difference between 
the two observed excitations is increased compared with 
what would have been observed for an isolated spin-1 
chain. We applied the correction for this Qa-dependence 



and the Qa-dependence 
q) to extract the disper- 



of the dispersion by using Eq. 
of the excitations (inset of Fig 
sion of an uncoupled chain. The 3D corrected dispersion 
is shown in the inset of Fig. ^ showing that the energy 
of the excitation towards Qc = is clearly higher than 
towards Qc = 1- 

The Qc-dependence of the excitation width (Fig. ||) 
suggests a cross-over from a one-particle dispersion for 
Qc > 0.5 to a two-particle continuum with the corre- 
sponding increased excitation width for Qc < 0.5. The 
cross-over from a one-particle to a two-paijticle response 
has also been observed by Zaliznyak et alSB at 1.5 K in 
the 3D ordered phase. In their experiment the increase 
in width is clearer because the thermal broadening of the 
one-particle peaks is much less. In both experiments the 
broadened peaks are centered at energies lower than the 
onset of the two-particle continuum for non-interacting 
chains (Fig. ||). Possibly this is because of the 3D inter- 
chain coupling which at 1.5 and 6.2 K gives considerable 
dispersion to the excitation energies for Qc = 1 as Qa 
varies. This could decrease the two-particle onset energy. 



subtracting phonon scattering near 10 meV, which was 
estimated from measurements at large wave- vectors, the 
integrated continuum intensity is 9(2)% of the total in- 
tensity at Q = (0.81, 0.|S 1,1), consistent with the previ- 
ous estimate of 12(2)%.eI 




Energy (meV) 

FIG. 7. Neutron scattering intensity at Q = (0.81,0.81, 1) 
(solid circles) as a function of energy transfer measured us- 
ing the DUALSPEC triple-axis spectrometer. The scattering 
intensity below 3 meV is shown 20 times reduced. The scat- 
tering at Q = (0.81,0.81,0.5) is shown for comparison (open 
circles). The solid and dashed line are fits to the measured 
spectra as explained in the text. 



B. Continuum scattering 

Fig. shows two constant- Q scans measured using the 
DUALSPEC triple-axis spectrometer at the Chalk River 
Laboratories and compares the scattering at the AF point 
Q = (0.81,0.81,1) with that at Q = (0.81,0.81,0.5), 
both at T = 6 K. For Qc = 1 new data with the 
sample mounted so that no hydrogenous material was 
in the beam are presented. The data are corrected for 
neutron absorption. This was done by measuring the 
quasi-elastic scattering when the sample presented the 
same angles to the incident and scattered beams as in 
the constant-Q scans. At Qc = 0.5, where the disper- 
sion of the well-defined excitation has a maximum, the 
scattering intensity below and above the sharp excitation 
is considerably lower than the intensity at Qc = 1, and 
it gives a good measure of the nonmagnetic background. 
The excess scattering at Qc = 1 at energies above the 
single-particle excitation is the continuum scattering. 

The well-defined excitation at Q = (0.81,0.81,1) 
is well described by an antisymmetrized Lorentzian 
weighted by the Bose factor and convoluted witlLlhe res- 
olution ellipsoid given by Cooper and Nathans,EZl but it 
does not account for the continuum scattering. After 



The wave- vector dependence of the multi-particle con- 
tinuum on Qc was measured at 6.2 K using the time-of- 
flight spectrometer MARI. Three constant-Qc scans at 
(Qc) = 0.7, (Qc) = 0.8 and (Qc) = 1 shown in Fig. | 
confirm that for energies higher than the energy of the 
well-defined excitation, a slowly decreasing intensity is 
observed for all three wave- vector transfers. 

The observed spectra were fitted by an antisym- 
metrized Voigt function weighted by the Bose factor and 
convoluted with the line-shape of the quasi-elastic inco- 
herent scattering peak, which should give a good account 
of the energy resolution at low energy transfers. A Voigt 
function, which is a convolution of a Lorentzian with a 
Gaussian function, is needed to account for the extra en- 
ergy width resulting from the wide wave- vector resolution 
AQc = 0.1. The fits give a good account of the well- 
defined excitations (Fig. ||), but do not account for the 
scattering at higher energies associated with the multi- 
particle continuum. For (Qc) — 1, the scattering extends 
up to about 12 meV. The scattering at (Qc) = 0.8 and 
(Qc) = 0.7 extends to slightly lower energy transfers, 
indicating a decreasing upper energy boundary of the 
continuum with increasing |Qc — 1|. 
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FIG. 8. Neutron scattering spectrum at 6.2 K for various 
wave-vector transfers. The peak intensities are scaled to fit 
onto the graph. The lower solid lines are fits to the peak data 
and the upper line is a guide to the eye enclosing the shaded 
area of the continuum scattering. The data were obtained on 
the MARI spectrometer with Ei — 30 meV. 

In order to confirm that the extra scattering is not a 
resolution artifact, a series of tests using a spectrometer 
simulation programme available at ISIS were made to in- 
vestigate vi^hethcr the continuum above the wcU-dcfincd 
excitations could result from the convolution of the neu- 
tron resolution function with the single mode dispersion. 
The program takes into account the spectrometer param- 
eters - the detailed pulse line-shape, the chopper charac- 
teristics, the position and dimensions of the detectors 
etc. - and predicts the scattering line-shape for a par- 
ticular sample orientation and model cross-section. The 
results of this simulation reproduced the measured line- 
shape of the quasi-elastic incoherent scattering and of 
the well-defined excitations, and showed that none of the 
scattering at higher energies could arise from the tail of 
the resolution function. 

We also investigated whether the scattering contin- 
uum can be influenced by neutron scattering by phonons. 
Measurements were performed at 8, 12 and 200 K, with 
an incident energy of Ei ~ 20 meV, and with the ID 
axis oriented perpendicular to the incident neutron beam 
direction. The scattering intensity, normalized to the in- 
cident beam monitor, is shown in Fig. |^ at different tem- 
peratures. 



FIG. 9. Upper panel: Neutron scattering spectrum at 
Qc = 1, at the average of T = 8 and 12 K (solid circles), 
and at T = 200 K (open circles). The intensity was inte- 
grated for 0.9 < Qc < 1.1. The data were measured using 
the MARI spectrometer and Ei = 20 meV. Lower panel: the 
imaginary part of the magnetic susceptibility x"(^,'^) as a 
function of energy transfer for the quantum-disordered phase 
at low temperatures and in the high-temperature limit. The 
grey solid circles give the estimated background as described 
in the text. 

Above 5 meV, the intensity of the continuum scatter- 
ing does not change between 6 and 12 K (Fig. |l^). We 
then averaged the spectra at 8 and 12 K to obtain better 
statistics at high energy transfers for (T) = 10 K. This 
leads to an artificial broadening of the well-defined exci- 
tations, and so the averaged data at low energy transfers 
were not analyzed in detail. 

The upper panel in Fig. ^ compares the neutron scat- 
tering intensity at Qc = 1 measured with MARI in the 
quantum disordered phase at (T) = 10 K (8 — 12 K) 
with that measured at 200 K. The lower panel of Fig. O 
shows the imaginary part of the generalized suscepti- 
bility x"iQc = 1,^^) = S{Qc = l,w)/(n(w) -1- 1) ob- 
tained from the dynamic structure factor. Yet at 200 K 
x"{Qc = l,w) is clearly less than at (T) = 10 K for all 
energy transfers, apart from a sharp peak at 11.5 meV 
which can be associated with a single-phonon excitation. 
It is not present in the measurements with incident en- 
ergy Ei = 30 meV which sampled a different scatter- 
ing surface. If the continuum scattering between 5 and 
12 meV for (T) = 10 K were due to neutron scatter- 
ing by phonons, x"{^^ ^) would be comparable at 200 K. 
The fact that it decreases greatly shows that the slowly 
decreasing scattering intensity above the Haldane mode 
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is magnetic and not phonon scattering. At T = 200 K, 
x"{Qc = 1,1^) has a broad peak at 'v-S meV, as we will 
discuss in a forthcoming publicationEl This broad peak 
is absent at low temperatures when only low-energy mag- 
netic excitations are thermally activated and x"(7r,a;) is 
dominated by the slightly-broadened Haldane gap mode. 
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FIG. 10. Neutron scattering intensity at T = 6.2 K as a 
function of wave- vector transfer along the chain, Qc, and in- 
tegrated between 8 and 12 meV (solid circles) and integrated 
between 12 and 16 meV (open circles). Lines are guides to the 
eye. The data were measured in all three detector banks of the 
MART spectrometer with an incident energy Ei = 30 meV. 



An estimate of the non-magnetic background was also 
obtained from the scattering observed in detectors at 
high scattering angles 2G for which the wave-vector 
transfer |Q| is large, and the magnetic form factor is 
low so that the spectrum is dominated by the phonon 
scattering which is proportional to |Qp or |Q|^.E3 The 
non-magnetic background for the low-angle scans was es- 
timated by scaling the high- angle intensity as |Qp. The 
background that is independent of |Q| was estimated 
from the energy gain side of the spectrum at low tempera- 
tures. The total background including phonon scattering 
is shown in Fig. ^ as grey solid circles. Assuming that the 
whole intensity at high scattering angles scales with |Qp 
or that part of the scattering scales with |Q|^ leads to 
a yet smaller background. The estimated non-magnetic 
background scattering accounts only for a small fraction 
of the observed scattering at Qc = 1 up to 12 meV and 
substantiates the magnetic origin of the observed contin- 
uum. 

The momentum dependence of the continuum shows 
a broad maximum at Qc — 1. Fig. ^ compares two 
constant-w scans for (w) = 10 meV and for (w) = 
14 meV; the latter can be taken as the background. The 
constant-tj data were obtained by integrating the inten- 
sity between 8 and 12 meV and between 12 and 16 meV, 
respectively. Thus the magnetic scattering from the well- 
defined excitation was not included in the integration. 
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11. Energy-integrated continuum scattering at 
T = 6.2 K as a fraction of the integrated intensity at 
(0.81,0.81,1) and as a function of wave-vector transfer Qc 
along the chain. The data is corrected for the magnetic form 
factor. The solid line is a guide to the eye. 

Fig. |l^ shows that for wave- vector transfers 0.6 < Qc < 
1.4, the scattering intensity for (w) = 10 meV is higher 
than for {lo) = 14 meV. The fact that the additional 
scattering is located around the AF zone center is further 
evidence for its magnetic origin. 

The energy-integrated intensity of the multi-particle 
states for 0.5 < Qc < 1 is shown in Fig. |l^ as a fraction 
of the total scattering intensity at Qc = 1 at the ID 
point (0.81,0.81, 1). It was determined numerically from 
the measured intensity after the integrated intensity of 
the well-defined excitation was subtracted. Fig. |ll| shows 
that the energy-integrated intensity of the continuum has 
its maximum at Qc = 1 and approaches zero for Qc = 
0.5. The integrated intensity is smaller near Qc = 0.5 
partly because the range of energies for which the multi- 
particle continuum can be observed gets smaller as the 
momentum decreases. 

Sum Rule: A further test of the magnetic na- 
ture of the continuum scattering is provided by the 
Hohenberg-Brinkman first moment sum rule. The first 
energy moment F{Qc) = J dujujS{Q,Lu) of ID magnets 
with nearest-neighbor exchange interaction is given by 
F{Qc) cx (1 — cos(7rQc)).cj This relation is valid evea. 
in the presence of a weak interchain coupling constant 
The first energy moment F{Qc) was determined from the 
observed spectra after subtracting a flat background and 
correcting for the magnetic form factor. The contribu- 
tion from the well-defined excitation is shown in Fig. |l^ 
as open circles and from the total scattering by the solid 
circles. The solid line is the prediction of the Hohenberg- 
Brinkman sum rule which predicts that the first moment 
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doubles between Qc = 0.5 and Qc = 1- It is clear that the 
observed spectrum is consistent with the sum rule, while 
F{Qc) determined from the well-defined excitations alone 
is not. This result confirms that the multi-particle con- 
tinuum is magnetic, that the single-mode approximation 
fails and that the presence of the continuum is required 
to fulfill the sum rule. 
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FIG. 12. i^(Qc) = J dLJujS{Q,uj) at T = 6.2 K as a func- 
tion of the wave-vector transfer Qc along the chain. Solid 
circles show F{Q) calculated numerically from the observed 
spectra, open circles represent F(Q) of the well-defined excita- 
tions. The solid line is the result of the Hohenberg-Brinkman 
sum rule. 



The presence of the high-energy continuum has impor- 
tant consequences for the static structure factor S{Qc)- 
Because scattering at higher energies enters the first mo- 
ment with a higher weight than low-energy scattering, 
S{Qc) is smaller if high-energy excitations are present. 
Since the continuum is most intense for Qc = 1 and lower 
for increasing |Qc— 1| the width of S{Qc) as a function of 
Qc is expected to increase around Qc = 1, corresponding 
to a reduced correlation length. This is what we have 
indeed observed (see later). 

Temperature Dependence: The temperature depen- 
dence of the multi-particle continuum was investigated 
between 6.2 and 12 K. Fig. 13 compares the scattering 
intensity at Qc = 1 for 6.2 and 12 K. The two spectra 
were measured in the same configuration, with the inci- 
dent energy Ei — 30 meV and the c-axis perpendicular 
to the incident neutron beam direction. At the higher 
temperature, the well-defined excitation has a higher en- 
ergy and becomes wider ,|--consistent with our previous 
triple-axis measurements. t3 The model of the antisym- 
metrized Lorentzian function weighted by the Bose fac- 
tor describes the well-defined excitations well, but the 
fits cannot account for the scattering at higher energies 
(Fig. At both T = 6.2 and 12 K the muhi-particle 



scattering extends up to about 12 meV. Between 6.2 and 
12 K the multi-particle scattering is largely independent 
of temperature and the scattered intensity between 5 and 
12 meV energy transfer is tlie same at both temperatures. 
This result was used earlierO to average the two data sets 
to produce a color plot with increased statistics. 
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FIG. 13. Neutron scattering intensity as a function of en- 
ergy transfer for wave-vector transfer with 0.9 < Qc < 1.1 
at 6.2 and 12 K. The peak from the Haldane gap mode is 
shown with its intensity reduced by a factor of 10. The lines 
are fits of a Lorentzian function to that observed intensity at 
6 K (dashed line) and 12 K (solid line). 



The upward renormalization of the excitation en- 
ergy and the broadening of the excitation with increas- 
ing temperature increases the scattering at higher en- 
ergy transfers as the temperature is increased. Much 
above 12 K the multi-particle states cannot be distin- 
guished unamhiffupusly from the broadened single parti- 
cle excitation.BEl 



C. Correlation length 

The instantaneous structure factor at T = 6.2 K for 
the ID chain, S{Qc), is shown in Fig. |lj. It was de- 
termined experimentally by integrating the scattering 
observed with the MARI spectrometer over both the 
well-defined excitations and the continuum scattering. 
The background scattering was subtracted, the remain- 
ing scattering corrected for the magnetic form factor de- 
pendence and then integrated numerically up to energy 
transfers of 12 meV to give S{Q). In the experiment the 
intensity was measured close to but not exactly at the 
non-interacting wave- vector (inset in Fig. The struc- 
ture factor of the independent chain was obtained from 
S{Q) by assuming as for antiferromagnets that the inte- 
grated intensity of the excitations close to the AF point 
is proportional to 1/lu{Q) and the ID structure factor is 
given by 
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uj{Qc) 



(7) 



where uj{Qc) and ti'(Q) are given by Eqs. || and ||. 

If the correlations along the chains decay exponentially, 
as for a ID Ising model, the structure factor is given by 



sinh(^^ 



cosh(^j ) + cos(7r(5c) 



(8) 



where the correlation length ^7 is measured in spin steps 
c/2 along the chain. If ^ ^ 1 and \Qc — 1| ^ 1 then this 
can be approximately written as a Lorentzian 



SiiQc) = 



26 



l + 7r2(g,-l)2^/ 



(9) 



When this function is convoluted with a Gaussian to take 
account of the wave- vector range sampled along the chain 
direction, and fitted to the S{Qc) shown in Fig. the 
correlation length obtained is 2.2(2) spins. It would be 
only 1.9(1) spins if the 3D correction, Eq. 0, is neglected 
because the data near Qc — 1 were taken for Qa < 0.19 
where S{Q) is suppressed. 

The theory of the structure factor of = 1 linear 
chains is usually formulated using the single-mode ap- 
proximation in which it is assumed that all the spectral 
weight resides in a mode whose energy coincides with 
the fast moment uj{Qc) of the Hohenberg-Brinkman sum 
ruleJHa 



4 <H>(l-cos^Q,) 



(10) 



where < Ti. > / L is the ground state energy per spin. Us- 
ing Eq. H for (jj(Qc), S{Qc) near Qc = 1 is proportional 
to a square-root Lorentzian 



'S'sm(Qc) oc 



where we expect 



A 



(11) 



(12) 



In Eq. |Tl| we neglect a small logarithmic term necessary 
to conserve the total spin. When our experimental values 
oivs, a and A are inserted we expect ^ — 4.8(2). 

The square-root Lorentzian gives a better description 
of the data in Fig, [l^ than the Lorentzian, Eq. ^ After 
convoluting Eq. gdj with a Gaussian to take account of the 
wave- vector range sampled for each Qc, and fitting to the 
data between 0.7 < Qc < 1-3, we find ^ = 4.0(2) sites and 
if the 3D correction, Eq. 0, is neglected 3.0(2) sites. The 
value of 4.0(2) less thait-f — 5 from numerical calcula- 
tions for isolated chainsH It is also less than the value of 
4.8(2) deduced from the parameters of the well-defined 
excitations indicating that the single-mode approxima- 
tion is not entirely valid. The correlation length deduced 



from Eq. ^ is shorter because of the presence of the con- 
tinuum, which depresses S{Qc) near Qc = 1 as may be 
seen from the first moment sum rule. Qualitatively the 
continuum excitations at high energy contribute less to 
S{Qc) than would a spectrum that was entirely confined 
to low-energy excitations. 
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FIG. 14. Static from factor, S{Qc), for the wave- vector 
transfer Qc along the chain as a function of energy trans- 
fer at T = 6.2 K. The open circles represent data measured 
at negative wave-vector transfer Qc and reflected at Qc = 0. 
The solid line is a fit to the data, based on Eq. |l^, as explained 
in the text. 



The observed value of ^ at 6.2 K is much lower than 
the result reported by Zaliznyak et alH This is be- 
cause their measurements were performed with wave- 
vectors close to the minimum of the excitation energy 
at Qoidcr = (0.33,0.33, 1). Consequently, the results are 
strongly influenced by 3D effects because the correlation 
length measured at this wave-vector diverges at the 3D 
ordering temperature T/v. Those results cannot therefore 
be compared with the correlation length of isolated quan- 
tum chains. Furthermore Zaliznyak et alB performed the 
energy integration only up to energy transfers of 1.6 meV 
and our results show that there is appreciable scattering 
at higher energy transfers. This effect and their chosen 
wave- vector transfer both give rise to a larger correlation 
length that that observed in our experiment. 



V. DISCUSSION 

With improved measurements of the gapped quantum- 
disordered spin system CsNiCla at 6.2 K we have mapped 
out the behavior of weakly-coupled spin-1 chains. We 
find (1) that the Haldane excitations have a spin veloc- 
ity Vs — 5.70(7) jaeV — 2.5 J - in excellent afiteement 
with field theoryl^j and numerical calculations,! 
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(2) that the single-chain correlation length is shorter than 
that of isolated chains, and (3) that the observed multi- 
particle continuum at the AF point carries about 12(2)% 
of the total scattering. This confirms previous experi- 
mental resultgj, that the continuum is much larger than 
expected theoretically. 

A. Continuum 

Multi-particle continua for AF wave- vectors have been 
predicted for T = K in the framework of the NLcrM and 
numerical techniques. In the framework of the NLtrM 
the continuum consists of three-particle states whose dy- 
namic structure factor has a energy gap of 3Ao and a 
pronounced maximum in its intensity at ^ GAn .(Ap is 
the Haldane gap energy at zero temperature). Its 
continuum integrated intensity from SAq up to 20Ao is 
only a tiny fraction, 1%, of the total scattering at Qc = 1, 
even if the coupling of the magnetic chains are takenJjito 
account in a Random Phase approximation (RPA)E3 

The diagonalization of the quantum eigcnstates of 
a spin-1 chain with = 20 sites predicts that at 
Qc = 1 the Haldane excitation carries 97% of the to- 
tal scatteringj23 leaving no more than 3% of the scatter- 
ing weight to contribute to continuum scattering. More 
precise results can be obtained using the density matrix 
renormalization group method. These calculations show 
that at T = K the Haldane excitation in a spin-1 chain 
with A^ = 256 sitp. carries 97.6(1)% of the total scat- 
tering at Qc = lj23 predicting thus at least 2.4 times 
more scattering than the NLcrM. This discrepancy is 
not very surprising because the three-particle states at 
higher energies involve excitations far away from the AF 
wave-vector where the field theory is known not to give 
accurate results. 

Clearly none of the theories can explain our ob- 
served result of 12% continuum scattering. We con- 
clude therefore that the multi-particle continuum ob- 
served in CsNiCla is not in accord with theoretical results 
for an isolated AF spin-1 Heisenberg chain with nearest- 
neighbor interactions at T = K. 

B. Comparison to Majorana Fermion Theory 

Other models for AF spin-1 chains assume that the 
spins also interact via biquadratic exchanges. Bi- 
quadratic spin interactions are thought to be small in 
most materials, but they cannot be ruled out a priori. 
The Hamiltonian for these models can be written as 

chain 

n = jY, Si.Si+i + fe(s,.s,+i)2. (13) 

i 

This Hamiltonian can be exactly solved for b — —1 (Ar- 
menian point) and for 6 = 1/3 which is the valence-bond- 



solid (VBS). It has been shown numerically that the en- 
ergy gap of this Hamiltonian i|S-|non-zero for — 1 < 6 < 1 
with a maximum at 6 = 0.41.E3 Close to the Armenian 
point, the spin system can be treateiin the framework 
of Tsvelik's Majorana fermion theorycj, which is a low- 
energy field theory for linear spin-1 chains. This theory 
predicts that the multi-particle scattering states jCarry 
17% of the intensity at Qc = 1 for b close to -l.Ea At 
the VBS point, the excitation spectrum has an energy 
gap A = 0.664J.tj Numerical calculations predict that 
the integrated intensity of the well-defined excitation at 
Qc — I carries 99% of the total scattering and multi- 
particle states are weaker by a factor of 2-compared to 
the Heisenberg point of the model (6 — 0).e2I 

In order to test whether the multi-particle continuum 
could be caused by a strong biquadratic term contained 
in the Hamiltonian of the Majorana theory, where 17% 
of the scattering at Qc — 1 lies in the continuum, we 
have searched at 6 K, where the gap is A = 1.24 meV, 
for the onset of the continuum at 3A and its pronounced 
maximumEj at ~ 5A. Our resolution is more than ade- 
quate to observe a dip between the one-particle peak at 
A and a maximum of the continuum at 5A = 6.2 meV, 
but instead the spectrum shows a steadily decreasing in- 
tensity in this energy range (Fig. |^, ^). This may arise 
because the 3D interactions cause a change in the en- 
ergy of the gap with the in-plane wave-vector Q^ At 
6.2 K, the minimum gap energy is about 0.4 mcVQ and 
so the minimum energy for the onset of the continuum 
is 1.2 meV and comparable with the ID gap energy A. 
This effect together with the thermal broadening of the 
one-particle excitations may be sufficient to fill in the gap 
between the one-particle and the continuum contribution 
to the scattering. 

More importantly, however, a strong biquadratic term 
(6=1) not only increases the spectral weight of the 
mukii-particle continuum, it also reduces the energy 
gapO in the excitation spectrum to zero. This is not in 
accordance with the facts. In CsNiCls the single-chain 
gap A at 6.2 K is higher than expected from the bilinear 
exchange J. This suggests that biquadratic interactions 
are weak and that the observed magnetic multi-particle 
continuum cannot be explained by the Majorana fermion 
theory. 

C. Correlation length 

The single-chain correlation length ^ has been deter- 
mined to be ^ = 4.0(2), which is considerably shorter 
than ^ = 5 obtained from quantum Monte Cpjlo calcula- 
tions for finite chains at T = 0.25J ~ 6.2 K.l^ Since the 
measurements of ^ were made close to the non-interacting 
wave-vector, Qa — 0.81, they should give the correlation 
length of isolated chains. Nevertheless the discrepancy 
may arise if the RPA approximation, Eq. |^, does not re- 
move all of the 3D effects adequately. 
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The ground state of an isolated chain has a hid- 
den string order. It is composed of states such as 

I H 1 — + 00 — 0+ > in which the antiferromagnetic 

ordering is maintained across intervening zero (Sz — 0) 
sites. Possibly this order is destroyed more effectively by 
3D interactions than is predicted by the RPA approxi- 
mation. 

The strong quantum fluctuations in = 1 linear chains 
are largely due to the low spin and low dimensionality. 
The fluctuations may further be enhanced, and the cor- 
relation length decreased, if there are frustrated inter- 
actions. In CsNiCla there are two possibilities for frus- 
trated interactions. Firstly the chains are arranged on a 
hexagonal lattice so there are bound to be frustrated in- 
teractions between the chains arranged at the corners of 
the triangle. A second possibility is that the inter-chain 
interactions may be between not only spins in the same 
basal plane but also between spins in the neighboring 
chains that are at ±c/2 along the chains.ciJ If these inter- 
actions are also strong and antiferromagnetic, as might 
be the case for super-exchange interactions, then even the 
interactions between pairs of chains are at least partially 
frustrated. 

Further theoretical work is needed to understand 
whether these multiply-frustrated interchain interactions 
are responsible for the shortened correlation length as 
well as tbe-associated upward renormalization of the Hal- 
dane gapcj, or whether the effects result from a failure 
of the RPA approximation. 

D. Origin of continuum 

We first compare with a material with smaller 3D in- 
teractions. Ni(C2H8N2)2N02C104 (NENP) is a AF spin- 
1 chain system in which the interchain couplings are not 
frustrated and much smaller than in CsNiCla , but a large 
single-site anisotropy splits tha.three Haldane excitations 
into two distinct excitations.c3 For NENP at 0.3 K the 
ratio of the first moment at the AF point and at the 
maximum of the dispersion is i? = F(l)/F(0.5) ^ 1.8(3) 
for the weighted average of the doublet^ and singlet ex- 
citations if data published by Ma et al.u is used for the 
calculationSy,and this is in agreement with the theoreti- 
cal value 2.c2l Neglecting the multi-particle continuum in 
CsNiCla this same calculation yields R = 1-2(2). This 
shows that in contrast to CsNiCls, the sharp excitations 
in NENP are enough to fulfill the Hohenberg-Brinkman 
sum rule and no magnetic weight is expected at higher 
energy transfers as in CsNiCla at 6.2 K. Any contin- 
uum would be hard to observe even in highly-deuterated 
NENP because of residual hydrogen scattering. Because 
NENP has weaker 3D interactions, this comparison sug- 
gests that interchain interactions may be responsible for 
generating the continuum observed in CsNiCla. 

We speculate that the multi-particle continuum is gen- 
erated by the AF coupling of the spin-1 chains, and mul- 



tiply frustrated interchain interactions might be decisive 
for the size of the effect. A shortened AF correlation 
length - due to an frustration-induced enhancement of 
quantum fluctuations - has two main effects on the dy- 
namic structure factor at Qc = 1. Firstly, it leads to 
an additional upward renormalization of the gap energy 
which is proportional to 1 (see Eq. |l^) . Secondly, a de- 
creased correlation length results in a reduced structure 
factor S{Qc ~ 1) because the total scattering must be 
constant (Eq. ^l]). Additionally there is strong evidence 
that the absolute value of S'(Q = (0.81,0.81, 1)) at low 
temperatures is smaller than that of isolated chains by at 
least 1.6£2l If only a well-defined mode existed it would 
have to increase its energy at Qc = 1 by a factor of 1.6 
to satisfy tlie first moment sum rule of Hohenberg and 
BrinkmanEa However, the Haldane energy is only 20% 
higher relative to an uncoupled chain and so the increase 
of Haldane energy is insufficient to compensate for the 
loss of intensity. The only way for the system to comply 
with the sum rule is to create a scattering continuum at 
high energies, as we find it does. This is evidence that 
the observed continuum in CsNiCla is correlated with an 
enhancement of quantum fluctuations due to multiply- 
frustrated 3D spin correlations. 

E. Similarity to two-spinon continuum of spin- 1/2 
chains 

Fig. ^ compares the observed intensity at 6.2 K (left 
panel) with the scattering expected, for an AF spin- 1/2 
chain based on the Miiller AnsatzCJ (right panel). The 
peak of the lower energy boundary of the two-spinon 
continuum was chosen to match the maximum of the 
single-particle dispersion in CsNiCla for direct compar- 
ison of the two excitation spectra. The experimentally 
observed integrated intensity /tot — J dQc J dto S{Qc, i^) 
of CsNiCla with S" = 1 was set equal to S{S -I- 1) = 2. 
The intensity of the Miiller Ansatz, which is valid for 
chains with S = 1/2, was scaled to 2S{S -I- 1) = 1.5 (the 
factor 2 accounts for the spin-1/2 components per Ni^+ 
site arising from the two holes in the d-shcU of the ion 

Ni2+). _ . . rn 

The two excitation spectra in Fig. ^ have obvious 
similarities: the multi-particle continuum extends to ap- 
proximately 12 meV, which is twice the maximum of 
the one-particle energy, and the upper boundary of the 
continuum decreases with increasing distance from the 
AF point. Both continuum spectra can be described 
with a lower boundary t^min — wafzbI sin((5c7'')| (the 
experiment could not establish whether the continuum 
scattering is gapped or not) and an upper boundary 
2wmax = t^AFZB sin(S^), whcrc WAFZB = 6 meV. Qual- 
itatively the continuum spectrum found in CsNiCla has 
similarj-features to the continuum found in AF S = 1/2 
chains. □ 
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FIG. 15. Left panel: Scattering intensity observed in 
CsNiCls at T = 6.2 K as a function of wave- vector trans- 
fer Qa and energy transfer on a logarithmic scale. The 
data were measured with an incident energy Ei = 30 meV. 
Data measured in detectors at high scattering angles are 
not shown. Right panel: Scattering intensity expected for 
a spin-1/2 chain according to the Miiller Ansatz as a func- 
tiop-iof wave-vector transfer and energy transfer (in units of 
J)m The energy axis of the Miiller Ansatz was scaled so that 
the lower bound of the two-spinon continuum matches the 
maximum of the single-particle dispersion in CsNiCla . 

The intensity of the continuum in a spin-1/2 chain is 
much stronger than that in CsNiCls. This is shown in 
Fig. |l^, which compares the scattering at Qc = 1 as a 
function of energy transfer for CsNiCla and for a spin- 
1/2 chain. The lines of equal intensity are also different 
in the two cases (Fig. |l5|). At 7 meV, as an example, 
the intensity in CsNiCls is constant, while it has a clear 
minimum at Qc = 1 for the Miiller Ansatz. 

We conclude that the intensity and the wave-vector 
dependence of the observed continuum suggests that a 
small amount of spin-1/2 degrees of freedom are gener- 
ated in CsNiCla, possibly by multiply- frustrated inter- 
chain interactions. In analogy to the excitations in a 
spin-1/2 chain, the excitations in CsNiCls may therefore 
also be described by a single spectral function containing 
a sharp onset and a continuum at higher energies. We. 
found that the spectral functions given by field theoryC3 
are indeed able to describe both the well-defined excita- 
tion and the multi-particle continuum close to the AF 
point if the critical exponents are adjusted, suggesting 
that the well-defined excitation is an onset of a scatter- 
ing continuum. More theoretical work is needed to clarify 
the spectral function and its dependence on the 3D inter- 
actions as the likely source of the multi-particle states. 




0.0 



2 4 6 8 10 12 14 
Energy (meV) 

FIG. 16. Measured scattering intensity of CsNiCla at 
(Qc) = 1 and at T = 6.2 K as a function of energy trans- 
fer. The solid line is the scattering intensity expected from a 
spin-1/2 chain according to the Miiller Ansatz. 



VI. CONCLUSIONS 

In summary, the excitation spectrum of CsNiCla has 
been measured between 6.2 and 12 K. The results con- 
firm that the spin system is in a quantum-disordered 
phase and the excitation spectrum has its full Qc = 2- 
periodicity. The single-chain correlation length at 6.2 K 
is 4.0(2) sites along the chain, which is lower than the 
value predicted by quantum Monte Carlo calculations, 
possibly due to spin correlations perpendicular to the 
chain axis. The reduced correlation length and the as- 
sociated upward renormalization of the Haldane gap un- 
derlines the need for more powerful theoretical tools to 
describe the frustrated AF coupling of quantum chains. 

The key result is the observation of a magnetic multi- 
particle continuum for wave- vectors 0.6 <Qc < 1.4 along 
the chain. The magnetic multi-particle states have higher 
energies than the well-defined excitations and extend up 
to about twice the maximum in the one-particle disper- 
sion. The integrated intensity of the multi-particle con- 
tinuum is much larger than that predicted for the three- 
particle continuum by field theoretical models and nu- 
merical techniques for AF spin-1 chains with nearest- 
neighbor exchange interaction. The observed multi- 
particle continuum resembles the two-spinon continuum 
of AF spin-1/2 chains, but its intensity is weaker. This 
suggests the presence of spin-1/2 degrees of freedom in 
the quantum-disordered phase just above the ordering 
temperature, possibly generated by competing interac- 
tions due to the 3D AF interactions. 
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